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We argue that a systems of tightly knotted, Hnked, or braided flux tubes will have a universal 
mass-energy spectrum, since the length of fixed radius flux tubes depend only on the topology of 
the configuration. We motivate the discussion with plasma physics examples, then concentrate on 
the model of giueballs as knotted QCD flux tubes. Other applications will also be discussed. 
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I. INTRODUCTION 

It is known from plasma physics that linked magnetic 
flux tubes are much more stable than an unknotted single 
loop Pi . Linked flux tubes carry topological charge, and 
this can be thought of as a conserved (at least to low- 
est order) physical quantum number. Similarly, knotted 
flux tubes carry topological quantum numbers, and one 
can think of a knot as a self-linked loop. The topological 
charges are described by knot polynomials that are re- 
lated to projections of knots or links into a plane where 
the crossings of the loops are assigned various attributes. 
Following each line around its loop generates the polyno- 
mials. Several types of polynomials have been studied in 
the literature (see e.g. Refs. 0,0): Alexander, Conway, 
Jones, Kauffman, etc., with increasing levels of precision 
for distinguishing knots. For example, the simplest knot, 
the trefoil, has a chiral partner (mirror image) that is 
not detected by the simpler polynomials, but is by the 
more sophisticated ones. Hence, a pair of knots with dif- 
ferent polynomials are different, but the converse is not 
necessarily true. It is still an unsolved problem to find a 
set of polynomials that distinguishes all non-isomorphic 
knots/links. Similar results hold for braids, and we will 
also discuss these objects below. 



II. REVIEW OF PREVIOUS PHYSICAL 
RESULTS ON TIGHT KNOTS AND LINKS 



If the loops have fixed uniform thickness and circu- 
lar cross-section (we will eventually discuss how one can 
relax this condition) , then each knot and link has a com- 
pletely specified length if the configuration is tight, i.e., 
is of the shortest length with the tubes undistorted and 
non-overlapping. If tubes have uniform cross sections, as 
can be approximately the case with magnetic or electric 
fiux tubes carrying quantized fiux, or for a polymer or 
even a piece of spaghetti, then the length of the tight 
knot is proportional to the mass (or energy) of the knot. 
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This, we claim, generates a universal mass (energy) spec- 
trum for knotted/linked configurations of objects of this 
type. The lengths of tight knots were not studied until 
the mid-1990s and only recently have accurate calcu- 
lations of large numbers of tight knots 5] and links [H 
become available. These results now make it possible to 
examine physical systems and compare them with the 
knot spectrum. The first physical example studied was 
tightly knotted DNA |^]. More recently, we have exam- 
ined the glueball spectrum of QCD . These particles Q 
are likely to be solitonic states '9'| that are solutions to the 
QCD field equations. While QCD will be our main focus 
in this chapter, there are many more cases where tight 
knots may play a role. We first proceed with an analysis 
of flux tubes in plasma physics. The lack of controllable 
quantum flux renders this case somewhat less interest- 
ing than its generalization to QCD. We will not go into 
any experimental details here, but we hope the experts 
in the areas discussed will take our general perspective 
into account when analyzing their data. 

In order to decide if a system of flux tubes falls into 
the universal class of having a tight knot energy spec- 
trum, we must flrst investigate the time scales involved. 
These are the lifetime of the soliton Ts and the relax- 
ation time Tr necessary to reach the ground state of a 
tight knot conflguration. The soliton lifetime (or the 
corresponding decay width F^ ~ ^/ts) can depend on 
several factors. These include the effects of flux tube 
breaking, rearrangement, and reconnection. The partial 
width for flux tube breaking is non-zero if the produc- 
tion of particle/anti-particle pairs is energetically possi- 
ble, for example monopole/anti-monopole (MM) pairs or 
color monopole/anti-monopole (McMc) pairs for mag- 
netic flux (or color magnetic flux) or quark/anti-quark 
(qq) pairs for color electric flux tubes. The partial widths 
can vary widely depending on the particle masses (e.g., 
iriq ^ M, so we expect qq pairs to be easier to produce 
than MM pairs), interaction strengths (this, for instance, 
enhances MM pair production versus qq pair produc- 
tion), and boundary conditions (tube shape and length). 
Rearrangement is a quantum effect where, for example, in 
a double donut arrangement, the loops can tunnel free of 
each other. Finally, reconnection is another effect where 
tubes break and re-attach in a different configuration. 
Such behavior has been seen in plasma physics, and is 
of major importance in understanding a variety of astro- 
physical systems. All these processes change topological 
charge, and their partial widths compete more or less 
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favorably with each other depending on the parameters 
that describe the system. 



III. EXACT CALCULATIONS 

While no knot lengths have been calculated exactly, it 
is possible to calculate the exact lengths of an infinite 
number of links and many braids . For links these cal- 
culations are possible in the case where each individual 
elements of the link lies in a plane. For braids, exact cal- 
culations are possible when the elements of the braid are 
either straight sections or where their centerlines follow 
helical paths. The shortest of all links, the double donut, 
is exactly calculable. The two elements lie in perpen- 
dicular planes and are tori of equal length. The shortest 
non-trivial braid is a helically twisted pair. "Weyl's tube 
formula" is the ideal tool for calculating the volume of 
a flux tube 0, 0| . The formula states that for a tube 
of constant cross-section a normal to a path of length /, 
the volume of the tube is just Vt = l<y- This is a remark- 
able result that holds in flat 3D. In higher dimensions or 
curved space, the result is more complicated, but here we 
need only the simple 3D case. This means that if we have 
an analytic form for the path and a circular cross-section 
we can find Vt- This leads to the fact that there is a class 
of exactly calculable links and braids. Since there are no 
known analytic forms for the path of tight knots, their 
volumes can only be calculated numerically, but once we 
have an estimate of the length, we can then also estimate 
the volume and therefore the energy for a corresponding 
physical system. Braids are also of physical interest. For 
the simplest example, a tightly twisted pair, the path 
of the center lines of such tubes are helices, and so the 
lengths l{9p,h) depend on the pitch angle 9p, the radius 
of the tubes r, and the height of the braid h. Hence the 
volume is Vt = 2al{6p, h) where a — irr^ . The volume of 
triple, etc., helical tight braids can also be found exactly; 
however, as with knots, the volumes of topologically non- 
trivial tight braids (those where the elements are woven 
together) can only be found approximately. (The helical 
twisted pair can also become non-trivially knotted/linked 
by identifying top and bottom boundaries, but we will 
not pursue this possibility here.) While the simple he- 
lically twisted braid has a volume that depends on the 
pitch angle which can potentially be adjusted by experi- 
mental conditions (see below), the tight knots and links 
derived from braids have no such adjustable parameter. 



IV. PLASMA PHYSICS 

Before going on to our main example of QCD, let us 
stop here to discuss tight links of flux in electromagnetic 
plasma. This example is conceptually somewhat simpler 
and provides motivation for what is to come. 

Movement of fluids often exhibits topological proper- 
ties (for a mathematical review see e.g. Ref. |l3)- For 



conductive fluids, interrelation between hydro- and mag- 
netic dynamics may cause magnetic fields, in their turn, 
to exhibit topological properties as well. For example, for 
a perfectly conducting fluid, the (abelian) magnetic he- 
licity J cPx e^^^ AidjAk is an invariant of the motion [l3l |. 
and this quantity can be interpreted in terms of knot- 
tedness of magnetic flux lines [Tjl . Let us discuss this in 
some detail since its implications are central to our more 
general results. 

A. Magnetic relaxation 

A perfectly conducting, viscous and incompressible 
fluid relaxes to a state of magnetic equilibrium without a 
change in topology P| . The system approaches a state of 
magnetic equilibrium by decreasing its magnetic energy, 
which is achieved by contraction of the magnetic field 
lines. In the case of trivial topology, where field lines are 
unknotted and unlinked closed curves that can be con- 
tracted to a point without crossing each other, such mag- 
netic relaxation proceeds uninterrupted. For example, 
the initial toroidal field configuration upon contraction 
deforms into a configuration of poloidal fields confined 
to a tube perpendicular to the original torus. Such con- 
figuration is still unstable since small disturbances aug- 
mented by the magnetic pressure lead to the increase of 
length of the tube and decrease of its cross-section. The 
relaxation eventually leads to a state with zero fields (vac- 
uum). 

If, however, the topology of the initial magnetic fields 
is non-trivial (for example, when flux tubes are knot- 
ted or linked), the relaxation stops when flux tubes are 
tightly knotted or linked. This happens because the 
"freeze-in" condition (see below) forces topological re- 
strictions on possible changes in field configurations and 
so any initial knots and links of field lines remain topo- 
logically unchanged during relaxation. The energy of a 
final (equilibrium) state is determined by topology and 
is bounded from below. One such bound is proportional 
to \H\, where TL is the magnetic helicity (see Ref. (ij for 
details). 

B. Abelian helicity 

Consider an abelian gauge potential 1-form A = 
Aidx^ and the corresponding field-strength 2-form F — 
^Fijdx^dx^ . The helicity for the field inside volume V is 
defined by 

[ AF. (1) 
Jv 

Under a gauge transformation, A ^ A + dip and F — > i^, 
and so using the Bianchi identity dF = and the Stokes 
theorem we find 

n^n + [ -ipF. (2) 
Jav 
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The helicity is thus gauge invariant if the normal compo- 
nent of the field F vanishes on the surface dV. 

It is easy to calculate the helicity for two linked flux 
tubes with fluxes $1 and <I>2. Considering first infinitely 
thin tubes (centered around the curves Ci and C2) and 
integrating over their cross sections wc find 

H = /" A + $2 /" A (3) 

The Stokes theorem now leads to 

n = 2n$i$2, (4) 

where n is the Gauss linking number of the two tubes 
(the algebraic number of times that one tube crosses the 
surface spanned by the other tube). It is straigtforward 
to generalize this to the case of linked and/or self-linked 
thick fiux tubes. 



C. Non-abelian helicity 

Now we begin to step toward QCD by considering a 
non-abelian plasma. Since both A and F are conserved 
during the plasma motion, any combination of A and F 
is a candidate for a conserved quantum number Ti.. The 
choice of expression is narrowed by requiring that 7i is a 
topological quantity. In particular, this requires that H 
be a surface integral. By analogy with the abelian case, 
for a conserved non-abelian helicity [l5|, we choose the 
corresponding expression with topological properties, 

[ {AdA + ^A"^). (5) 
Jv 

A (time-independent) gauge transformation A — > 
g~^Ag + g~^dg leads to 

n-^n-^tif {g'^dgf + tif gdg-^A. (6) 

Jv JdV 

We consider only a limited set of gauge transformations 
such that do not change the condition v^Af^ = 0; this im- 
plies the constraint v^dig — 0. Because of this constraint 
only two of the components g~^dig are independent and 
so (g^^dg)^ = 0. The helicity is invariant if the normal 
component of A vanishes on the surface dV. This condi- 
tion is more restrictive than the one needed in the abelian 
case and implies the latter. 

D. "Freeze-in" condition 

In a perfectly conductive relativistic non-abelian 
plasma, the electric field vanishes in the local frame mov- 
ing with the plasma. The Lorentz transformation to the 
rest frame gives v'^ F^u = 0, where = (1, w*) is the 
local plasma velocity and 

i^M- = - d^A^ + [A^, A,] (7) 



is the field strength in terms of the gauge potential A^ . 
(The equation v^F^^i, = is the appropriate general- 
ization of its familiar abelian non-relativistic counter- 
part ^5 Ei = CijkV^B^.) In the "hydrodynamic gauge" 
v^Af^ = 0, from the definition of the field strength F^^, 
we obtain 

doA, + {dtA.y + Akid^v") = (8) 

and 

doF,, + {OkF^^y + F,k{d,v^) + i^fe,(a,t.") = 0. (9) 

As a fluid particle moves in space, the rate of change 
of a local quantity is given by the Lagrange derivative 
d/dt ~ v^di_i.. Applying this operator to differential forms 
A = Aidx^ and F = ^Fi-jdx^dx^ , and using Eqs. ® 
and 

{d/dt)dx' = {d.jv^)dx^, (10) 

we find that A and F are constants of motion: 

dA/dt = 0, (11) 
dF/dt = 0. (12) 

In other words, for a curve C and a surface S moving 
with the fiuid, integrals /^^ A and F do not change in 
time: the magnetic lines are "frozen" into the plasma. 
Note the integral A is usually called a Wilson loop in 
the high-energy physics literature and F is related to 
the Wilson loop via Stokes' theorem. 

The freeze-in condition derived in this subsection gen- 
eralizes the well-known (see e.g. Ref. jlQ]) "freeze- 
in" condition in MHD to non-abelian case. The non- 
relativistic analog of this result was obtained in Ref. ■ 

V. QCD AND GLUEBALLS 
A. QCD 

What is the ideal physical system in which to discover 
and study tight knots and links? We claim it is Quan- 
tum Chromodynamics (QCD), so in order to explain our 
reasoning wc first pause to briefly summarize QCD Q. 

To lowest order, the standard model of particle physics 
can be broken into two major sectors, one describing 
the electroweak processes for leptons, quarks, photons 
and intermediate vector bosons, like electron-neutrino 
scattering, and the other describing the strong interac- 
tions. The electroweak processes are easier to deal with 
because the coupling constants of this sector are small 
(e.g., the fine structure constant a = e^/47r is approxi- 
mately 1/137) so that perturbative calculations can be 
preformed to a high precision and are in beautiful agree- 
ment with experiment. On the other hand, the strong 
interactions which describes the interactions of quarks 
and gluons, while in principle completely described by 
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QCD, are much more difficult to deal with because at 
low energy the couphng constant as = gs^/47r is 0(1). 
This, for instance, makes the quark-antiquark {qq) bound 
state problem analytically intractable. 

In somewhat more detail, the QCD is a gauge theory 
with the lagrangian density 

^ - -i^;.^'^'^"^ + i^\i^{D^),,i^i - m,r,r,, (13) 

where the field strength is 

f;, - d^Al + d.A'l ~ gsfabcAlAl (14) 
and the covariant derivative is 

(D^^j = %<9^ + (15) 

Here the tp^^s are the Dirac spinor quark fields with 
color (i) and flavor (q) indices, the Yang-Mills fields 
A° (connection) describe the gluons, the fate are the 
structure constants for the SU{3) gauge group, and the 
A° 's are the triplet representation matrix elements for 
the quarks. Quarks come in six flavors with six associ- 
ated quantum numbers conserved by the strong inter- 
actions but all except electric charge can be violated 
by the electroweak interactions. Color charges are con- 
fined due to the fact that SU (3) is an unbroken symme- 
try with asymptotic freedom. I.e., the strong coupling 
as is energy dependent and becomes small at high en- 
ergy but is large at low energies where it becomes 0(1) 
at a few hundred MeV. Then the theory becomes non- 
perturbative and confining (color charges cannot be iso- 
lated) with all observables color singlet states. These 
singlets are described by the quark model and are either 
bosonic qq bound state called mesons (pions, kaons, etc.) 
or fermionic qqq states called baryons (protons, neutrons, 
etc.). Besides the mesons and baryons there are a number 
states that do not fit neatly into the quark model. These 
go by the names of hybrid states, exotic states and glue- 
balls. Hybrid states and exotic states are thought to be 
unusual combinations of quarks e.g., qqqq bosonic states, 
or qqqqq fermionic states, most of which have net flavor 
charges. The glueballs on the other hand are thought 
to be made of gluons with at most some virtual quark 
content, hence no flavor charges. 

We are now in a position to support our claim that 
QCD is the ideal physical system in which to discover 
and study tight knots and links. Here are our reasons: 

1. QCD is a solidly based part of the standard model 
of particle physics, and much about color confine- 
ment and the quark model is already well under- 
stood in this context, making much previous work 
transferable to the problem of tightly knotted flux 
tubes in this theory. 

2. Unlike plasmas, fluids or other condensed matter 
systems where flux tubes are excitations of some 
media with many parameters that could hide uni- 
versal behavior, flux knots in QCD can exist in the 



vacuum. Thus continuum states are absent and 
there are no media parameters to vary and obscure 
the universality. Hence, the results in QCD can be 
far less ambiguous. 

3. The hadronic energy spectrum has been measured 
over a large range of energies (140 MeV to 10 GeV) 
and already many hundreds of states are known. 
We expect that among these, a few dozen can 
be classified as tightly knotted/linked flux tubes 
states. These states must have no valance quarks 
(i.e., no flavor quantum numbers) in order to be 
classified as glueballs. 

4. Knotted solitons in QFT are already known to ex- 
ist. 

5. One can efficiently search for new glueball states at 
accelerators. (Also, data from older experiments 
still exist and can be reanalyzed to check the pre- 
dictions of new states described below.) 

B. Knot energies 

Consider a hadronic collision that produces some num- 
ber of baryons and mesons plus a gluonic state in the 
form of a closed QCD flux tube (or a set of tubes). From 
an initial state, the fields in the flux tubes quickly re- 
lax to an equilibrium conflguration, which is topologi- 
cally equivalent to the initial state. (We assume topolog- 
ical quantum numbers are conserved during this rapid 
process.) The relaxation proceeds through minimization 
of the fleld energy. Flux conservation and energy min- 
imization force the fields to be homogeneous across the 
tube cross sections. This process occurs via shrinking 
the tube length, and halts to form a "tight" knot or 
link. The radial scale will be set by Aq^j-,. The en- 
ergy of the final state depends only on the topology of 
the initial state and can be estimated as follows. An ar- 
bitrarily knotted tube of radius a and length I has the 
volume na'^l. Using conservation of flux ^e, the energy 
becomes oc Z(tr<I>|;)/(7ra^). Fixing the radius of the tube 
(to be proportional to AqJ,j-,), we find that the energy 
is proportional to the length I. The dimensionless ratio 
e{K) ~ I /{2a) is a topological invariant and the simplest 
definition of the "knot energy" [l8| . 

Many knot energies have been calculated by Monte 
Carlo methods and certain types can be calculated 
exactly (see below) , while for other cases simple estimates 
can be made (see TablelJ). For example, the knot energy 
of the connected product of two knots Ki and K2 satisfies 

e{Ki#K2) <e{Ki)+e{K2). (16) 
A rule of thumb is 

s{Ki#K2) ~ e{Ki) + e{K2) - (2^ - 4). (17) 
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TABLE I: Comparison between the glueball mass spectrum and knot energies. 



oXiSiXiQ 


iviass 


i\ 






/o(600) 


400 - 1200 


2! 


12.6 [47r] 


768 [766] 


/o(980) 


980 ± 10 


3i 


16.4 


993 


/2(1270) 


1275.4 ± 1.2 


2? * Oi 


[67r + 2] 


[1256] 


/i(1285) 


1281.9 ± 0.6 


4i 


21.2 


1277 






A 1 

4! 


(21.4) 


(1289) 


/i(1420) 


1426.3 ± 1.1 


5i 


24.2 


1454 


{/2(1430) 


« 1430} 


5i 


24.2 


1454 + 
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1200 - 1500 


3i * Oi 


(24.7) 


(1484) 
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1507 ± 5 


52 


24.9 


1496 
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1518 ± 5} 


52 


24.9 


1496 + (5 


/a (1525) 


1525 ± 5 


52 


24.9 


1496 + 35 


{/2(1565) 


1546 ± 12} 


r 2 

51 


(25.9) 


(1555) 


{/2(1640) 


1638 ± 6} 


63 


((27.3)) 


((1638)) 






(2\ *0i) *0i " 


[87r + 3] 


[1686] ^ 


/o(1710) 


1713 ±6 


6i 


((28.6)) 
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3i#3| 


28.9 (30.5) 
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3i#3i 


29.1 (30.5) 


1744 (1827) 






2f * 2{ 


[87r + 4] 


[1745] 






62 


29.2 


1750 






61 


29.3 


1756 






63 


30.5 


1827 






7i 


30.9 


1850 






819 


31.0 


1856 






820 


32.7 
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/2(2010) 


201li^E! 


72 


33.2 
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/4(2050) 


2025 ± 8 


821 


33.9 


2028 






81 


37.0 


2211 






10l61,162 


37.6 


2247 


/2(2300) 


2297 ± 28 


818, 9l 


38.3 


2288 


/2(2340) 


2339 ± 60 


92 


40.0 


2389 






lOi 


44.8 


2672 






111 


47.0 


2802 



"Notation nj^ means a link of I components with n crossings, and 
occurring in the standard table of links (see e.g. 0) on the fc**^ 
place. K^K' stands for the knot product (connected sum) of knots 
K and K' and K * K' is the link of the knots K and K' . 

''Values are from |^ except for our exact calculations of 2\, 2\*0\, 
and (22 *0i)*0i in square brackets, our analytic estimates given in 
parentheses, and our rough estimates given in double parentheses. 

=E(G) is obtained from e(K) using the fit TT^ . 

"^States in braces are not in the Particle Data Group (PDG) sum- 
mary tables. 

■^This is the link product that is not 2^ * 2j. 

■^Resonances have been seen in this region, but are uncon- 
firmed 0. 

Most of the knot energies in Table 1 have been taken 
from Ref. but we have independently calculated the 
energy of 2\, A\ and 6| exactly and the energy for several 
other knots and links approximately. We find £(2^) — 



Att « 12.57, to be compared with the Monte Carlo value 
12.6. We also find e(4?) = Gtt -I- 2 and e(6f ) = Stt -)- 3, 
where there are no Monte Carlo comparisons available, 
or needed. Other exactly calculable links can be found 
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in Ref. 6"! and an example of a link with energy IOtt + 5 
is shown in Fig. ^ 

C. Model 

In our model, the chromoelectric fields Foi are confined 
to knotted/linked tubes. After an initial time evolution, 
the system reaches a static equilibrium state which is 
described by the energy density 

£e = ^tvFo.F''' - V. (18) 

Similar to the MIT bag model we have included a 
constant potential energy V needed to keep the tubes 
at a fixed cross-section. The chromoelectric flux is 
conserved and we assume flux tubes carry one flux quan- 
tum. To account for conservation of the flux, we add to 
Eq. lO the term 

tvX{^E/i7ra^)-n'Fo,}, (19) 

where is the normal vector to a section of the tube and 
A is a Lagrange multiplier. The energy density should 
be constant under variations of the degrees of freedom, 
the gauge potentials A^. This leads to the system of 
equations 

i?"(Fo, -An,) =0, (20) 
D'{Fo^-Xn,) = 0, (21) 

which has the constant field 

Fo^ - ($B/™2)n, (22) 

as its solution. With this solution, the energy is positive 
and proportional to I and thus the minimum of the energy 
is achieved by shortening I, i.e., tightening the knot. 

The case of chromomagnetic flux tubes can be simi- 
larly considered. This requires the confinement of color 
magnetic flux tubes which is possible if there are light 
quarks in the spectrum of the theory [2(1 |. 

Lattice calculations, QCD sum rules, electric flux tube 
models, and constituent glue models agree that the light- 
est non-qq states are glueballs with quantum numbers 
J++ = 0++ and 2++ (see Ref. IH). We wiU model 
all 7++ states (i.e., all fj and f'j states listed by the 
PDG 0), some of which will be identifled with rotational 
excitations, as knotted/linked chromoelectric QCD flux 
tubes. We proceed to identify knotted and linked QCD 
flux tubes with glueballs, where we include all /,/ and 
fj states. The lightest candidate is the /o(600), which 
we identify with the shortest knot/link, i.e., the 2^ link; 
the /o(980) is identified with the next shortest knot, the 
3i trefoil knot, and so forth. All knot and link ener- 
gies have been calculated for states with energies less 
then 1680 MeV. Above 1680 MeV the number of knots 
and links grows rapidly, and few of their energies have 
been calculated. However, we do find knot energies cor- 
responding to all known fj and f'j states, and so can 



make preliminary identifications in this region. (We fo- 
cus on fj and f'j states from the PDG summary tables. 
The experimental errors are also quoted from the PDG. 
There are a number of additional states reported in the 
extended tables, but some of this data is either conflicting 
or inconclusive.) 

Our detailed results are collected in Tabled where we 
list fj and f'j masses, our identifications of these states 
with knots and the corresponding knot energies. 

In Fig. 121 we compare the mass spectrum of fj states 
with the identified knot and link energies. Since errors for 
the knot energies in Ref. 0| were not reported, we con- 
servatively assumed the error to be 1%. A least squares 
fit to the most reliable data (below 1680 MeV) gives 

E{G) = (23.4 ± 46.1) + (59.1 ± 2.1)£(A') [MeV], (23) 

with = 9.1. The data used in this fit is the first seven 
fj states (filled circles in Fig. ^ in the PDG summary 
tables. Inclusion of the remaining seven (non-excitation) 
states (unfilled circles in Fig.[2Jl in Tabled where either 
the glueball or knot energies are less reliable, does not 
significantly alter the fit and leads to 

£;(G) = (26.9 ±24.9) + (58.9 ±1.0)e(A:) [MeV], (24) 

with = 10.1. The fit (|23|l is in good agreement with 
our model, where E{G) is proportional to s{K). Better 
HEP data and the calculation of more knot energies will 
provide further tests of the model and improve the high 
mass identification. 

For comparison, wc have fitted the same data in Fig.|21 
except that the first glueball is missed out; this results 
in = 231. Similarly, we have fitted the same data in 
Fig. ^except that this time the first link is missed out; 
this results in x"^ = 355. This is strong evidence that our 
identification is appropriate. 

In terms of the bag model the interior of tight 
knots correspond to the interior of the bag. The flux 
through the knot is supported by current sheets on the 
bag boundary (surface of the tube). Knot complexity 
can be reduced (or increased) by unknotting (knotting) 
operations Q, Q . In terms of flux tubes, these moves 
are equivalent to reconnection events [23] • Hence, a 
metastable glueball may decay via reconnection. Once 
all topological charge is lost, metastability is lost, and 
the decay proceeds to completion. Two other glueball 
decay processes are: flux tube (string) breaking; this fa- 
vors large decay widths for conflgurations with long flux 
tube components; and quantum fluctuations that unlink 
flux tubes; this would tend to broaden states with short 
flux tube components. As yet we are not able to go be- 
yond providing a phenomenological fit to these qualita- 
tive observations [7| , but hope to be able to do so in the 
future. 

We have assumed one fiuxoid per tube. There may be 
states with more than one fluxoid, but these would pre- 
sumably have somewhat fatter flux tubes with higher flux 
densities and higher energies. For example, the two flux- 
oid trefoil knot 3i would certainly have e{K) > 2e(3i) 



FIG. 1: An example of an exactly calculable link. 




FIG. 2: Relationship between the glueball spectrum E[G) and knot energies e{K). Each point in this figure represents a 
glueball identified with a knot or link. The straight line is our model and is drawn for the fit 1231 . 
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FIG. 3: The same as in Fig. ^ but tlie first glueball is missed out 
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and a fairly reliable estimate gives e{K) w 2V2e{3i). 
Hence most multifluxoid states would be above the mass 
range of known glueball candidates. 



VI. DISCUSSION AND CONCLUSIONS 

We expect all tight knots and links to be described 
by smooth curves with minimum radius of curvature 2r 
where r is the tube radius. Consider how we would ap- 
proximate such a curve on a square lattice. The simplest 
nontrivial example is the double donut. We consider one 
of its loops (a circle) which we call T to be centered on a 
particular lattice point P. The radius of the circle R is 
some number of lattice spacings away from P (we assume 
we are in a lattice plane), along the direction of one of 
P's nearest neighbors. Assume all other points on F are 
at a distance d > R away from P. If we minimize the 
length L of F, then the ratio of L to the circumference of 
a circle of radius R is 



which is the best approximation (about a 27% error) we 
can have on a square lattice. Similar arguments show we 
can do somewhat better on a hexagonal lattice where we 
can achieve 



27ri? 



= 2V3/TT, 



(26) 



2ttR 



A/tt, 



(25) 



about a 10% error. However, neither of these approx- 
imations are stunningly successful, nor does decreasing 
the lattice spacing improve the approximations. We as- 
sume other knots and links will be approximated on the 
lattice with similar accuracy. It is interesting to note 
here that both lattice estimates are too high and that 
lattice QCD typically predicts glueball masses above 1 
GeV, while Monte Carlo predictions for tight knots leads 
to a lighter glueball spectrum with lightest state just be- 
low 800 MeV. Lattice QCD are certainly more sophisti- 
cated and involved than the simple estimates we have just 
made, but we do expect requiring knots/links to lie on a 
lattice to give glueball masses on the high side. Further- 
more, the tube cross section on a lattice is not circular, 
inhibiting tight packing, and the amount of curvature en- 
ergy on the lattice is likely to be more than for a smooth 
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curve. We expect these effect to also contribute to in- 
creasing the lattice glueball mass estimates. 

Let us return to continuum physics and consider a slab 
of material that can support flux tubes. We have in 
mind a superfluid or superconductor, but are not lim- 
ited to these possibilities. Assume further that the flux 
tubes carry one and only one unit of flux. Next consider 
manipulating these flux tubes. For instance, consider 
a hypothetical superconductor where the flux tubes are 
pinned at the bottom of the slab, say by being attracted 
to the poles of some magnetic material, and at the top of 
the slab they are each associated with the pole of a mov- 
able permanent magnet, perhaps a magnetic whisker, or 
fine solenoid. If the north poles of a flux tube is on the 
top surface, we call it a flux tube (and (— ) for the 
south pole at the top). Since the flux tubes are pinned 
on the bottom of the slab we can maneuver the top ends 
into a braid. For example, if there are two flux tubes we 
can rotate them around each other to make a helically 
twisted pair. If there are three or more flux tubes we can 
twist them or braid them. If we measure the forces on 
the magnets when the flux tubes are being manipulated 
and use these to find the energy needed to form a braided 



configuration, then we gain information on the tension of 
the flux tubes and on the energy stored in a tight braid. 
If we have both and (— ) flux tubes we can bring 
- (— ) pairs together, annihilate the ends so the tube can 
be pulled into the interior of the superconductor. If we 
do this at both the top and bottom surfaces for a braided 
pair then it should be possible to form a bulk knot (sim- 
ilar for links). Assuming they have time to relax to tight 
configurations, the energy released by the eventual de- 
cay of these structures in the bulk should correspond to 
the universal energy spectrum described above. Perhaps 
there is no system with all the properties of our ideal 
superconductor; however, there are systems that contain 
many if not most of the desired features. 

Another collection of physical systems of potential in- 
terest for its ability to support vortices, and so knots and 
braids, are the atomic Bose-Einstein condensates. For 
example, dilute ^^Rb atoms at 80 nK. Recently results 
have been presented demonstrating laser stirring of these 
condensates to produce vortices [23. As many as seven 
vortices have been seen in a 0.1 mm^ region. One could 
hope that advancement in these techniques could lead 
to more complicated solitonic structures, in particular 
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knots, links, and braids. The advantage of these systems 
is the experimental accessibility of these structures, and 
the large percentage of the volume of the system taken up 
by the solitons. This could potentially dramatically im- 
prove the signal-to-noise ratio in the study of tight knots 
in condensed matter systems. 

We have argued that knotted/linked solitonic physi- 
cal systems have universal mass-energy spectra and have 



demonstrated this with a detailed example from QCD. 
We have provided other examples that are good candi- 
dates but where the universality has yet to be seen. A 
system is a candidate if it contains line solitons that can 
(somewhat paradoxically) relax to their tight knot/link 
ground state in a time shorter than their decay time. Fu- 
ture work along these lines could tie knot theory to many 
subfields of physics or to other sciences. 
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